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Abstract 

Up today estimations of proton-proton total cross sections, crj'^j, at very high energies are obtained 
from cosmic rays (> 10" eV) by means of some approximations and the knowledge of the measured 
proton-air cross section at these energies. Besides, crfj^ are measured with present day high energy 
colliders up to nearly 2 TeV in the center of mass (~ 10^^ eV in the laboratory) . Several theoretical, 
empirical or semi-empirical parameterizations, very successful for interpolation at accelerator energies, 
can then be used to extrapolate the measured value to cosmic ray energies and get a reasonable 
estimation of cross sections at higher energies (~ 10^^ eV). When both cross section estimations - 
from accelerators and cosmic rays - are compared, a disagreement is observed, amounting to more 
than 10%, showing a discrepancy beyond statistical errors. Here we use a phenomenological model 
based on the Multiple-Diffraction approach to successfully describe data at accelerator energies. Then 
we estimate with this model a^^^ at cosmic ray energies. The employed fitting method associates the 
model free-parameters to the data of only two physical observable: the differential cross-section and 
the parameter p. The model predictions for estimations of aj'^j are then compared with data of total 
cross-sections. On the basis of regression analysis we determine confident error bands, analyzing the 
sensitivity of our predictions to the employed data for extrapolation. : using data at 546 and 1.8 
TeV, our extrapolations for crfj^ are only compatible with the Akeno cosmic ray data, predicting a 
slower rise with energy than other cosmic ray results and other extrapolation methods. We discuss our 
results within the context of constraints in the light of future accelerator and cosmic ray experimental 
results. 



INTRODUCTION 



Recently a number of difficulties in uniting accelerator and cosmic ray values of proton-proton, cfoj , and 
antiproton-proton, crfot, hadronic total cross-sections, within the frame of the highest up-to-date data have 
been summarized P . Such united picture appears to be of the highest importance for the interpretation 
of results of new cosmic ray experiments as the HiRes and in designing proposals that are currently 
in progress as the Auger Observatory j^, as well as in designing detectors for future accelerators, as the 
CERN pp Large Hadron Collider (LHC) 0|. Although most of accelerator measurements of crf^j and 
crf^( at center of mass energies -/s < 1.8 TeV are quite consistent among them, this is unfortunately 
not the case for cosmic ray experiments at ^/s > 6 TeV where some disagreements exist among different 
experiments. This is also the case among different predictions from the extrapolation of accelerator data 
up to cosmic ray energies: whereas some works predict smaller values of crf^j than those of cosmic ray 
experiments (e.g. other predictions agree at some specific energies with cosmic ray results (e.g. |7|). 

Dispersion of cosmic ray results are mainly associated to the strong model-dependence of the relation 
between the basic hadron-hadron cross-section and the hadronic cross-section in air. The later determines 
the attenuation length of hadrons in the atmosphere, which is usually measured in different ways, and 
depends strongly on the rate (fc) of energy dissipation of the primary proton into the electromagnetic 
shower observed in the experiment: such a cascade is simulated by different Monte Carlo technics implying 
additional discrepancies between different experiments. Furthermore, ct^qj in cosmic ray experiments is 
determined from Up'^fair using a nucleon-nucleon scattering amplitude which is frequently in disagreement 
with most of accelerator data . 

On the other hand, we dispose of many parameterizations (purely theoretically, empirical or semi- 
empirical based) that fit pretty well the accelerator data. Most of them agree that at the energy of the 
future LHC (14 TeV in the center of mass) or higher the rise in energy of crf^j will continue, though the 
predicted values differ from model to model. We claim that both the cosmic ray and parametrization ap- 
proaches must complement each other in order to draw the best description of the proton-proton hadronic 
cross-section behavior at ultra high energies. However, the present status is that due to the fact that 
interpolation of accelerator data is nicely obtained with most of parametrization models, it is expected 
that their extrapolation to higher energies be highly confident: as a matter of fact, parameterizations 
are usually based in a short number of fundamental parameters, in contrast with the difficulties found 
in deriving crf^t from cosmic ray results Pj. With the aim of contributing to the understanding of this 
problem, in this paper we first briefly analyze in the first two sections the way estimations are done for 
cr^Qj from accelerators as well as from cosmic rays. We find serious discrepancies among both estima- 
tions. In section HI we describe the Multiple Diffraction model and the employed method to evaluate 
the model parameters on basis of data of only two physical observable: the differential cross-section and 
the parameter p. For the goal of the present approach we neglect crossing symmetry. In the fourth 
section we present a suitable parametrization to high energies of the free energy-dependent parameters 
of the model and we discuss its physical significance. In the fifth section we describe the procedure for 
the determination of error bands on basis of Appendix A. In section VI, on the basis of the Multiple 
Diffraction model applied to accelerator data, we predict a^^^ values with highly confident errors. We 
discuss in section VII our results in terms of the hypothesis crf^j — crf^j , and finally, we conclude with a 
discussion about the implications of extrapolations within the frame of present cosmic ray estimations. 
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1 HADRONIC aZ FROM ACCELERATORS 

Since the first results of tlie Intersecting Storage Rings(ISR) at CERN arrived in the 70s, it is a weU 
established fact that af^^ rises with energy (13 El ) ■ The CERN SppS Collider found this rising valid for 
(jfot as well 123- Later, the Tevatron at Fermilab confirmed that for crf^j the rising still continues at 1.8 
TeV, even if there is a disagreement among the different experiments as for the exact value 1121 [Qj . 
A full discussion on these problems may be found in |14l I15| . It remains now to estimate the amount 
of rising of the total cross section at those energies. Let us resume the standard technique used by 
accelerator experimentalists [H]. 

Using a semi-empirical parametrization based on Regge theory and asymptotic theorems experimen- 
talists have successively described their data from the ISR to the SppS energies. It takes into account 
all the available data for af^^, pPP, crf^j and ff^, where pPP'PP^ is the ratio of the real to the imaginary 
part of the {pp,pp) forward elastic amplitude at t = 0. The fits are performed using the once-subtracted 
dispersion relations: 

p^{EME)=- + -j^ (1) 

where Cg is the substraction constant. The expression for dtot'^ is: 

ai°V = A.E-''' ± A^E-^' + Co + C2Ms/so)]^ (2) 

where - (+) stands for pp (jjp) scattering. Cross sections are measured in mb and energy in GeV, E being 
the energy measured in the lab frame. The scale factor so has been arbitrarily chosen equal to 1 GeV^. 
The most interesting piece is the one controlling the high-energy behavior, given by a ln?{s) term, in 
order to be compatible, asymptotically, with the Froissart-Martin bound |16j . 

The parametrization assumes afj^^ and crfoj to be the same asymptotically. This is justified from the 
very precise measurement of the ppp parameter at 546 GeV at the SppS collider, ppp = 0.135 ±0.015 |17| . 
which implies that at present there is no sizeable contribution of the odd under crossing part of the forward 
amplitude, the so-called "Odderon hypothesis". This hypothesis predicts a value of ppp > 0.17— 0.20 

CHI, uni. 

The eight free parameters are determined by a fit which minimizes the function 

(7 pp Ppp f^pp Ppp ^ ' 

The fit has proved its validity predicting from the ISR pp and ppdata (ranging from 23 to 63 GeV in the 
center of mass), the af^f. value [201 later found at the SppS Collider, one order of magnitude higher in 
energy (546 GeV) JOI ■ With the same well-known method and using the most recent results it is possible 
to get estimations for af^^ at the LHC and higher energies. These estimations, together with our present 
experimental knowledge for both crf^j and crf^^ are plotted in figure 1. We have also plotted the cosmic 
ray experimental data from AKENO (now AGASSA) |2I1 and the Fly's Eye experiment [221 123 • The 
curve is the result of the fit described in [H] . The increase in crf^j as the energy increases is clearly seen. 

Numerical predictions from this analysis are given in Table 1 . It should be remarked that at the LHC 
energies and beyond the predicted crf^j and af^^ values from the fitting results display relatively high error 
values {^<yp°red — 8m6). Also, recent extrapolations for calculations of error bands, based on analysis 
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y = 2.2 (best fit) 
+-1 o 
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Figure 1: dg^'p from accelerator and from cosmic rays: solid line indicates the best fit obtained from a 
semi-empirical parametrization [H]. The two dashed lines delimitate the region of uncertainty. [ Cosmic 
ray estimations after \'21\ . 



Table 1: af^^ data from high energy accelerators: fits values are from [Hj. 



yCs (TeV) 




'y'Z (nib) 


0.55 


Fit 


61.8±0.7 




UA4 


62.2 ± 1.5 




CDF 


61.5± 1.0 


1.8 


Fit 


76.5 ±2.3 




E710 


72.8 ±3.1 




CDF 


80.3 ±2.3 




E811 


71.7±2.0 


16.0 


Fit 


111. 0± 8.0 


40.0 


Fit 


130.0 ± 13.0 
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taking into account cosmic ray and accelerator data, found error values relatively high, as for instance 5 
- 7 mb at 14 TeV [S]. We conclude that it is necessary to look for ways to reduce the errors and make 
the extrapolations more precise. 



2 HADRONIC aZ FROM COSMIC RAYS 

Cosmic ray experiments give us crf^j in an indirect way: we have to derive it from cosmic ray extensive air 
shower (EAS) data. But, as summarized in pp and widely discussed in the literature, the determination 
of crfoj is a rather complicated process with at least two well differentiated steps. In the first place, the 
primary interaction involved in EAS is proton-air; what it is determined through EAS is the p-inelastic 
cross section, crf^^"^, through some measure of the attenuation of the rate of showers. A™, deep in the 
atmosphere: 

A — i.\ — 14.5TOp 

— K^p—air — k p—air ^ ' 

The k factor parameterizes the rate at which the energy of the primary proton is dissipated into elec- 
tromagnetic energy. A simulation with a full representation of the hadronic interactions in the cascade 
is needed to calculate it. This is done by means of Monte Carlo simulations [251 126L 177) . Secondly, the 
connection between cr^^""" and afoj is model dependent. A theory for nuclei interactions must be used. 
Usually is Glauber's theory |28[ I29| . The whole procedure makes hard to get a general agreed value for 
crf^j. Depending on the particular assumptions made the values may oscillate by large amounts, from as 
low as 122 ± 11 at = 30 TeV quoted by the Fly's Eye group [22 ESI or 133 ± 10 mb by the Akeno 
Collaboration j^, both at = 30, to 162 ± 38 mb nearly = 30 ( around 160 - 170 mb ^/s = 40 
TeV), 1301 and even as high as 175l27also at yfs = 40 TeV |3J . In the last cases, and even taking into 
account the large quoted errors, the values for a^^^ are hardly compatible with the values obtained from 
the extrapolations of current accelerator data. 

From this analysis the conclusion is that cosmic ray estimations of crf^j are not of much help to 
constrain extrapolations from accelerator energies. Conversely we could ask if a trustable extrapolation 
based on accelerator data could not be used to constrain cosmic ray estimations. 



3 A MULTIPLE-DIFFRACTION APPROACH FOR EVALU- 
ATION OF aZ 



Let us tackle the mismatching between accelerator and cosmic ray estimations using the multiple- 
diffraction model applied to hadron-hadron scattering j32| . The elastic hadronic scattering amplitude 
for the collision of two hadrons A and B, neglecting spin, is described as 



F{q,s) 



I bdb 


1 _ e^m-s) 


Jo 





Jo{qb) 



(5) 



where ^(6, s) is the eikonal, h the impact parameter, Jo the zero-order Bessel function and — —t the 
four-momentum transfer squared. In the first Multiple Diffraction Theory the eikonal in the transferred 
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momentum space is proportional to the product of the hadronic form factors Ga and Gb (geometry) 
and the averaged elementary scattering amplitude among the constituent partons / (dynamics), and 
can be expressed at first order as ^{b,s) = Ca,b {GaGbJ), where the proportional factor Ca,b is the 
free parameter known as the absorption factor and the brackets denote the symmetrical two-dimensional 
Fourier transform. A connection between theory and experiment may be obtained by means of hadronic 
factors and elementary parton-parton amplitudes, which are not physical observable. However with the 
help of the optical theorem a^^f may be evaluated in terms of the elastic amplitude F{q, s): 



aZ^^^ ImF{q = 0,s) (6) 

crfot is a physical observable, the other two physical observable are the differential cross section and p 
expressed respectively as 

da 

and 



^^2_^|F(g,.)p (7) 



^ ReF{q^Q,s) 
^ ImF{q^O,s) ^ ' 

Multiple-diffraction models differ one from another by the particular choice of parameterizations made 
for Ga and Gb and the elementary amplitude /. In the case of identical particles, as is our case, 
Ga = Gb = G. For our purposes we follow the parametrization developed in a highly outstanding work 
[7] which has the advantage of using a small set of free parameters, five in total, which are in principle 
energy dependent: two of them {a^,0^) associated with the form factor 

G-(l + 4)"'(l + C)"' (9) 
and the other three (C, a. A) associated with the elementary complex amplitude / 

/(g,s) = i?e/(g,s)+*/m/(g,s) (10) 



where 
and 

so that, 



/m/(g,,s) = C^4j (11) 



Refiq,s) = Xis)Imfiq,s) (12) 
ImF{q 0, s) = 

1 - e-"('''") cos {Xn{b, s)} b db Jo{q, b) |,=o (13) 
with the opacity 0(6, s) given as: 

/•oc 

n{b, s)^ G" Imfiq, s) Jo{q, b) qdq (14) 

J O 
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which cxphcit expression is: 

nib, s) = C{EiKo{ab) + E^KoiPb) + E^K,,{ab)+ 

EiKeriab) + b [E^Ki{ab) + E^Kiipb)]} (15) 

where kQ,ki,kei, and ker are the modified Bessel functions, and Ei to Eq are functions of the five free 
parameters. The proton-proton total cross-section is directly determined by the expression 

aZ^'^TT J^bdb{l-e-''^'''''>cos[Xn{b,s)]'jjo{q,b) |,=o (16) 

expression that was numerically solved in |34| . It should be noted that, according to the principle 
of Analyticity, the scattering forward amplitude for particle-particle and particle-antiparticle appears 
from the same analytical function . The behavior of total cross-sections of both reactions is assumed 
to follow one of the following ways: up to the ISR energies the differences are attributed to Regge 
contributions, which are expected to disappear at higher energies ■ or on the contrary, differences 

are interpreted in terms of the Maximal Odderon Hypothesis, which predicts that they increase as the 
energy overpass the highest energy of the ISR. At this level let us now emphasize an essential point 
in all which follows. First, we quote the success in the prediction of the value of aZ at the SppS 
Collider JOI from the ISR data (mainly pp) using expressions (1) and (2), where aZ and aZ were 
taken asymptotically equal EOj. Secondly, according to |^ the difference Act between and CTf^j, 
Act = CT^Qj — CT^Qj, up to energies < 2000 GeV in the laboratory (~ 60 GeV in the center of mass) tends 
toward zero as s~° '^®. Figure 2. And finally, even if it is argued that and are anyway different 
for higher energies, we have indicated in Section I that current evidence points the other way round. 
Taking this triple line of evidence into account, in our multiple-diffraction analysis it is assumed the same 
behavior for aZ and aZ at high energy. So, here after = = ftot- It is noteworthy that some 
parametrization models, as the RRP ^3 predict the same value for both cross sections at y/s > 70 GeV 
and the same value for the corresponding p at y/s > 110 GeV. 

3.1 Evaluation of the Model Parameters 

For the evaluation of the free parameters it was assumed, following [Jj that two of them are constants: 
= 8.2 GeV^ and = 1.8 GeV^. The other three parameters are: C(s) and a~^(s) which determine 
the imaginary part of the hadronic amplitude [Eq. (11)] and consequently the total cross section [Eq. 
(6)] , and on the other hand, A(s), which controls the real part of the amplitude [Eq. (12)]. Setting A(s) 
— makes the amplitude purely imaginary, so that a zero (a minimum) is produced in the dip region, 
where only the real part of the amplitude becomes important We fit experimental data (Fig. 3) of 
the elastic differential cross section, with Eq. (7), choosing the set (C, a~^) for which the theoretical 
and the experimental central values are equated at the precise —q^ value where the data show its first 
minimum (the "dip" position), which coincides with the minimum of the imaginary part of the elastic 
amplitude. Data of pp differential cross section at 13.8 and 19.4 GeV was taken from ^U], for 23.5 — 62.5 
GeV from [HI, for 546 GeV (pp) from g^, (for - i < 0.5 (GeV/c)^) and [Sj (for 0.5 < ~t < 1.55 
(GeVC)2), and for 1800 GeV from |ll](for 0.034 < ~t < 0.65 (GeVC)^). The procedure is carried out for 
each one of the energies where there is available accelerator data of ^ in the interval 13.8 < ^/s < 1800 
GeV, as illustrated in Fig 4 and Fig 5. Because error bars of data at the dip position are small the 
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Figure 2: Plot of Aatot = o'p'p — Cpp* ^ s " ''^^ indicating that A(T*°* tends towards zero up to energies ^ 
60 GeV in the center of mass. Data points from [57] 



fitting procedure is based on the central values. It must be emphasized that, it is precisely because the 
minimum of the imaginary amplitude is produced in the dip region, that data at 1800 GeV is quite 
suitable to our procedure, since the predicted minimum falls a.t t = 0.585 GeV^ where there is available 
data = 0.0101). The next data value at t = 0.627 GeV'^ shows an slight increase. Furthermore, the 
shift of the dip region toward lower values of t as the considered energy increases is qualitatively consistent 
with the expectation relative to the shoulder at 546 GeV {t = 0.9 GeV^). Once the values of C(s) and 
a~^{s) are determined for each energy they are introduced, together with the constant parameters 
and into Eq. (8). Giving values to the parameter A(s) up to obtain the experimental value of p(s), 
we determine the A value for each energy where there is available accelerator data in the same interval 
13.8 < y/s < 1800 GeV. The obtained central values of the three energy-dependent free parameters, C(s), 
a~^(s)and A(s) are listed in Table 2. Therefore, the employed method to evaluate the model parameters 
only requires data of and p{s). 

4 PARAMETRIZATION FOR EXTRAPOLATION TO HIGH 
ENERGIES 

Concerning the physics behind a parametrization and fitting procedures, it should be emphasized that 
there is not yet an answer on microscopic basis, i.e., on terms of the QCD theory, because diffractive 
physics is essentially a non-perturvative phenomenon and confinement is the unsolved problem of QCD. 
However, on basis to the present model it can be said that the closeness of our fits to the colliders data 
(ISR, SppS and Tevatron) gives us an indication that the basic tenets of the model are right, as was 
described for instance in reference [32] . In spite of the approximation made in ignoring crossing symmetry. 
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Table 2: Values of the parameters C, a ^, A at each energy. They are obtained by equating the accelerator 
data and the model prediction for the elastic differential cross-sections and for the parameter p in the 
interval 13 < < 62.5 and 546 < ^/i < 1800 GeV. 





C{s) 


a-2(s) 




(GeV) 


(GeV-2) 


(GeV-2) 




13.8 


9.97 


2.092 


-0.126 


19.4 


10.05 


2.128 


-0.043 


23.5 


10.25 


2.174 


0.025 


30.7 


10.37 


2.222 


0.053 


44.7 


10.89 


2.299 


0.079 


52.8 


11.15 


2.350 


0.099 


62.5 


11.42 


2.380 


0.115 


546 


16.90 


3.060 


0.182 


1800 


21.52 


3.570 


0.194 




■ 1 1 2 3 4 5 e 7 a 9-1011 



|tj GeV' 

Figure 3: Data of pp differential cross section at ^/s < 62.5 according to 0ni,03i ESi 02] and 01]. 
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Itl(GeV') 

Figure 4: Fits of the predicted imaginary amplitude to data of the pp differential cross sections to 
determine C(s) and a~'^{s) by equating experimental and theoretical values at the specific | t | of the 
"dip" position for ^ = 546GeV and ^/s = ISOOGeV. 




j 12 3 4 6 6 7^ i. :,f Ifi 



Figure 5: Idem that Fig. 4, for ^/s = 52.8 and ^/s = 62.5GeV. 
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the parametrization carried out here is not an arbitrary one as wiU be described below. 

On this basis, with the aim of interpolating values at energies where there is no available accelerator 
data and for extrapolation to higher energies, we proceeded to perform a parametrization of the three 
energy-dependent free parameters. According to the obtained values of those parameters, as described in 
Section III-A (Table 2) and following the procedure to be described in Section V, a second-order fitting 
of the values of C(s) and a^^{s) and a exponential fitting of A(s) has been obtained from the following 
analytical expressions: 

C(s) = 19.24521 - 2.86114 ln{s) + 0.22616 ln{sf (17) 
a'^{s) = 1.8956 - 0.03937 Ini-s) + 0.01301 ln{sf (18) 

A(s) = 0.01686 + 0.00125 (l - e-;«(s/4oo)/o.i8549^ ^ 

0.19775(1- e-'"(V400)/3.74642^ (^g) 

Results are displayed in Table 2 and illustrated in Figures (6) -(8) as the central solid lines. 
However, it must be taken into account that the reliability of the functional parameterizations for extrap- 
olations must have a physical support, since it is clear that several parameterizations that may describe 
correctly the set of data in the experimental range 13.8 < ^/s < 1800 GeV, may not be consistent any 
more, but differ substantially when extrapolated to high energies. Therefore, some restrictions must 
be imposed in the selection of parameterizations according to the physical information available. As 
we mentioned before, our fittings of C{s) and a~'^{s) in the limited experimental range were based on 
experimental data of the differential cross-section and the p{s) parameter, obtaining values that increase 
with energy (Table 2) with positive concavity (Figs. 6-7). Experimentally, total cross-sections increase 
with energy like ln(s) or ln?{s) in the concerned energy range, and on the other hand, soft processes 
are expected to have a ln(s) behavior: from the optical theorem, the interdependence between the free 
parameters and the physical observable (Eq. 6) may be connected with the unitary condition, for which 
lowest order cross-sections within the frame of gauge field theories present ln(s) terms 45,. Based on 
our fittings and extrapolation results and on the previously mentioned constraints, the presence of ln(s) 
terms in the two-energy dependent parameters C(s) and Q!~^(s) is a natural consequence, and therefore 
the hypothesis of polynomial functions of ln(s) seems a reasonable one. 

Concerning the parameter A(s) let us pointed out that a basic property of the Glauber Multiple 
Diffraction Model is to associate elastic scattering cross-sections of nucleons with the scattering amplitude 
of their composite partons 32 : within this framework, the proportionality between real and imaginary 
parts of the parton-parton amplitudes (Eq. 12), A(s) — Ref{q,s)/Imf{q,s) implies that A(s) at the 
partonic scale behaves similarly as p{s) does at the hadronic scale. The influence of A(s) on the hadronic 
amplitude has been empirically analyzed in [7], showing that if A(s) increases (or decreases) also p{s) 
increases (or decreases) (Fig. 9), and A(s) = at the same energy value where p{s) = 0. However, due 
to the lack of p{s) data above the experimental energy range used in this work, the parametrization of 
A(s) at high energies is based on the conventional assumption that beyond y/s ^ 100 GeV, p{s) presents 
a maximum and then goes asymptotically to zero |14| . the rate of convergence depending on model 
particularities. Based on this consideration and on the empirical behavior of A(s), shown in Table 2, we 
propose the parametrization given in Eq. (19), where So = 400 GeV^ is the value at which A(s) converges 
to zero, and the numerical coefficients control the maximum and asymptotic behavior. 



10 



Let us now remind that blackening and expansion are very well known properties of elastic scat- 
tering. Within the context of the present empirical analysis, blackening and expansion are related to 
the elementary parton-parton amplitude and the hadronic form factors through the energy-dependent 
parameters C{s) and a{s) respectively. Since in the very forward direction the scattering amplitude 
is basically of diffractive nature, the eikonal becomes a purely imaginary one, then for hadron-hadron 
^(6, s) = C{s) {G'^Imf{q, s)) = Imil,{b,s), so that in terms of Eqs. (9) and (11) the opacity satisfies 
^{b, s) > 1, and the free parameter C(s) behaves as an absorption factor (optical theorem). On the other 
hand, the free parameter a{s)'^ may be connected to the hadronic radius [7] as i?^(s) = —6[dG{q, s)/dt]t=o. 
Then, from Eq. (9) 

1 1 



R^{s) = 0.2332 



a2(s) |3^s) 



ifmY 



Therefore, from Eq. (18) and the adopted value /S^ — 1.8 GeV^ the radius is an increasing function 
of energy and such a behavior is connected to the expansion effect. The result is that hadrons become 
blacker and larger with energy increase, consistently with the so called "Bell Effect" ^Hl- Since the 
hadronic scattering amplitude is purely imaginary the free parameters may be associated with the physical 
observable by means of Eq. (6) - (7). 



5 THE EXTRAPOLATION PROCEDURE 

The procedure followed to arrive to the predictions of crf^j at high energies with statistic confidence 
intervals based on the forecasting method described in Appendix A is the next: -(i) using the values 
displayed in Table 2, we established prediction equations of the type of eq. (A4) within the data range 
and eq. (A6) ) out of the range, for each of the energy-dependent parameters: (ii) using the least squares 
method in matrix formalism, as described in the Appendix (eqs. A10-A12) we obtained the constants 
00, $1, $2 for each one of the parameters. The autocorrelation constant was determined as described 
in ^||3S|, (iii) with the obtained values the central values of the three parameters were derived (eqs. 
17-19), leading to a second-order fitting of the values of C{s) and a^'^(s) and a exponential fitting of 
A(s); results are shown in Table 2, (iv) we evaluated the variance for each one of the new forecasted 
value (eq. A14) , (v) by means of eqs. (A15)-(A16) we estimated the confidence intervals for each one 
of the interpolated-extrapolated values, such that by a fitting of the extreme values of these confidence 
intervals we have built the error bands of each one of the parameters, what is shown in Figs. (6)- (8), 
(vi) the Central values of crf^j for each point are obtained by means of the introduction in eqs. (15)-(16) 
of the values of C(s), a~'^{s) and A(s) (displayed in Table 2), (vii) finally, the overall confidence band 
for the predicted trf^t is obtained, not from eqs. (A15)-(A16), but from the substitution of the extreme 
values of the error bands of the three energy-dependent parameters into eqs. (15)-(16), followed by the 
corresponding fittings (Figs. 10-11). 



6 RESULTS 

Data of the total cross sections with their respective errors are summarized in Table IV. At y/s < 62.5 
GeV data was taken from 41 , for 546 GeV from ^U]. For the value at 1800 GeV there exist three 
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Figure 6: The parameter C(s) from Table 2 with its confidence interval. 



Table 3: Fitted values (interpolation and extrapolation) of C(s), a ^(s) and A(s) . 





C{s) 
(GeV-2) 






(GeV) 


(GeV-2) 




13.8 


9.9039 


2.0945 


-0.12816 


19.4 


10.082 


2.1469 


-0.02848 


23.5 


10.225 


2.1798 
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Figure 9: Data behavior of the parameter p(s) at ^/s < 62.5 



different measurements: the value of the CDF collaboration (80.3 ± 2.24m6) |(12j, the value of the E710 
collaboration (72.8 ± 3.1mb) ^T] and the value of the E811 collaboration (71.7 ± 2.02m6) There is 
a high controversy on whether the correct value is that of the CDF collaboration or those lower values 
of the E710 collaboration and the E811 collaboration, mainly related to the estimation of the different 
backgrounds. For instance the CDF Collaboration has decided to use only their a^^i value to quote 
luminosity values for all of their physics program (SUI whereas the other Tevatron Collaboration, DO, has 
adopted the average of the three measurements For the goal of this work, one of us (J-V.) suggested 
to follow this second option and use the arithmetic weighted mean of the three values, (74.91 ± 1.35 mb). 
The consequence in our results of the employment of the values of those three collaborations is mentioned 
in the next section. 

We quote in Table IV the values of the predicted a^^^ for two different energy intervals: crj^t'^, 
(13.8 < < 1800 GeV), and crf^e, (13.8 < < 546 GeV). Figure 10 represents graphically the first 
case together with cosmic ray data. For further discussion, we include in Table IV predicted values from 
two standard accelerator-based data extrapolations: the first one, 0-540 makes a fit to all af^j and af^^ 
data in the interval 10 < y/s < 546 GeV using expression (2) and the second one, crlf' (Fig. 1), in 
the same energy range, makes the simultaneous fit to all crfot, ffot, Ptot ^^'^ Ptot data through the method 
described in Section I, . 
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Table 4: Central, upper and lower values for atot obtained with the forecasting method. The o'l^^'^, ^tot^ 
crf^j^ columns includes data up to 1800, 546, 62.5 GeV respectively. The quoted cr54g and crjf' values are 
from and|Hl respectively. 



7i (GeV) g-J (mb) aj^ (mb) 4^^' (mb) af^^ (mb) af^ (mb) ag,^ (mb) 

13^8 38.36 ± 0.04 38. 29^^;^^ ~ 38.3Qtoit ~ 38.5ll^ j^ 

19.4 38.97 ± 0.04 SSmto fr - - - 38.82tj] - - - SSMtoil 

23.5 38.94 ± 0.17 39.44+°^^ - - - 39.44+°j^ - - - 39.26toil 
30.7 40.14 ± 0.17 40.34l:o:95 - - - 40.37i^;|^ - - - AO.lbtom 

44.7 41.79 ± 0.16 Aimtl fe " " " 42.00t^|^ - - - 41.93lE;;^^ 

52.8 42.67 ± 0.19 42.76lo:J^^ - - - 42.84t°:^| - - - 42.93l^J^ 
62.5 43.32 ± 0.23 43.67+J;gg - - - 43.77+0:97 - - - 44.05lo 72 (*) 
546 61.5 ± 1.5 61.62+J:^^ - - - 61. 78+^:^^ - - - 69.39+^;4 

1800 74.91 ±1.35(*) 76.17ti:E;7 76.7 ± 4.0 76.00^2:34 76.5 ± 2.3 91.741^1? 

14000 - - - 108.2713^7 112 ± 13.0 106.5t^:^^ " " " 143.86t^^:^ 

16000 - - - 110.671^:™ - - - 108. 7t^:^:^ 111.0 ± 8.0 147.85t;^°:? 

30000 - - - 122.41^^:^^ - - - 119.6t|^^ - - - 167.64^:^^:^ 

40000 - - - 128.051^:°? - - - 124. 7t^:^? 130.0 ± 13.0 177.23t^^:J 

100000 - - - 147.141?:^^ - - - 142.0tii^:f - - - 210.06t^^:*J 

This value is the weighted arithmetic mean of the E710 (72.8 ± 3.1 mb), CDF (80.3 ± 2.3 mb) and E811 
(71.7 ± 2.0) values. 



7 DISCUSSION 



To start with, let first examine what happens when our main assumption, the asymptotic equality of 
crfoj and crfot, is not used. Analysis of Figure 11 shows that if we limit our fitting calculations to the 
accelerator domain ^/s < 62.5 GeV, where af^^ data exist, the extrapolation at cosmic ray energies 
produces an error band so large that practically any cosmic ray result become compatible with results 
at accelerator energies. It can be seen that in this case the af^^ values obtained when extrapolated to 
ultra high energies seem to confirm the highest quoted values of the cosmic ray experiments |3UI I31| . 
Also it can be noted that such extrapolation to ultra high energies may claim not only agreement with 
the analysis carried out in jSD] and the experimental data of the Fly's Eye |3T] , but even with the Akeno 
collaboration |21| , because their experimental errors are so big that they overlap with the errors reported 
in E|j and of course faUing within the predicted error band for that case ^/s < 62.5 GeV (Fig. 11). 
If that were true, it would imply the extrapolations cherished by experimentalists are wrong. But the 
prediction shown in Fig. 11 gives crf^j = 69 mb at the CERN SppS Collider (546 GeV), and 91.6m6 at 
the Fermilab Tevatron (1.8 TeV). Comparing with Table I we see that the measured af^j at 546 GeV is 
smaller than the predicted af^^ by near 8 mb, and in the case of 1.8 TeV by more than 15 mb, which no 
available model is able to explain |15j . However when our main hypothesis, af^^ — crf^j asymptotically 
is used, then the existing af^^ data at higher accelerator energies may safely be included. This allows 
to enlarge by a great amount the lever arm for the extrapolation and both the predicted values and the 
error band change considerably. This can be clearly seen in Fig. 10, as well as in Table IV, where we have 
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Figure 10: Prediction of the total cross-section using data up to 1800 GeV (solid line), (jlot'^- The error 
band (dashed lines) is obtained by the forecasting method. 



added the available trf^j up to 0.546 TeV (the af^^ column) and up to 1.8 TeV (the al^^° column). Now 
the predicted value of a^^^ from our extrapolation a^^^^, for instance at y^s — 40 TeV, (rf^j = 128.015 27 
mb, is seen incompatible with those in |H0[ K^lj by several standard deviations, though no so different to 
the Fly's Eye or Akeno results and the predicted value in 0. 

Concerning the quoted error bands, in sections V, it has been shown in the Appendix that due to 
the inclusion of Residual Correlations in the Forecasting method, this is always a high precise technique 
whatever the employed parametrization model, but this systematic advantage occurs, provided we are 
dealing with the same parametrization model under the same input conditions. Therefore, it must be 
emphasized that it is not valid a comparison between different statistical techniques on basis to different 
parametrization models, or running of the same parametrization with different input values. However, in 
spite of this, we would like to make a consideration of strictly qualitative nature, in the sense that other 
parameterizations 'SP with similar input quantities (for instance p and ^) lead to central values 
that are only slightly higher than ours, whereas the quoted errors are larger than ours by nearly a factor 
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Figure 11: Prediction of the total cross-section using pp data up to 62.5 GeV (solid line). The error band 
(dashed lines) is obtained by the forecasting method. 



of 3, as can be seen in Table 4, or in Fig. 4 and Fig. 10. For instance, the quoted error in the fit al^jr" at 
100 TeV in the center of mass (which corresponds to ~ 10^^ eV in the lab), llT^y gg, is comparable (or 
even better) to the error obtained at much more lower energies in other works as can be seen in Table IV, 
at 16 TeV and 14 TeV (~ 10^'' eV in the lab), corresponding to the energy range of future LHC CERN 
Collider. Another parametrization which predictions of the cental values are not so different of ours, 
is that of the Ensemble A, of the work Obviously, the previous analysis is only of qualitative 

nature, but it is however quite important because such a discussion opens new interrogatives: - though, 
it is clear that from statistic theory the forecasting technique is a highly precise method (section V, 
Appendix and references therein) , its advantage is, however, only a second order effect when applied to 
the same conditions (as shown in Eqs. 20 and 27), why this higher preciseness seems to be amplified when 
qualitative comparisons are made with a different parametrization model the same different inputs?, is 
it because the employed parametrization model is better?. It is clear that the answer must be based 
on a deep quantitative analysis, by applying the forecasting method to different parametrization models 
of total cross-section, in order to discern which model is more confident, and from there to draw more 
reliable physical inferences. Such an analysis is out of the scope of the present paper, and will be the 
subject of a forthcoming work. 
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8 CONCLUSIONS 



It has been shown in this work that highly confident predictions of high energy a^^^ values are strongly 
dependent on the energy range covered by experimental data and the available number of those data 
values. In particular, if we limit our study of determining af^j at cosmic ray energies from extrapolation 
of crf^j accelerator data in the range ^/s < 62.5 GeV, then the obtained results are compatible with most 
of cosmic ray experiments and other prediction models, because the predicted error band is so wide that 
covers their corresponding quoted errors (Fig. 11). However, as the included data in our calculations 
extends to higher energies, including af^^ data up to 1.8 TeV, that is, when all experimental available data 
is taken into account, the range of possibilities is much reduced. It should be noted that our predictions 
are compatibles with other prediction studies ^6,. Taken all this results at face value, i.e. as indicating 
the most probable af^^ value, we conclude that if predictions from accelerator data are correct, hence 
it should be of great help to normalize the corresponding values from cosmic ray experiments, as for 
instance by keeping the (k) parameter as a free one. The k value found will greatly help the tuning 
of the complicated Monte Carlo calculations used to evaluate the development of the showers induced 
by cosmic rays in the upper atmosphere. If extrapolations from our parametrization model are correct 
this would imply that o'pl'^^j^ should be smaller than usually considered, which would have important 
consequences for development of high energy cascades. Though it is quite clear that the statements of 
this work are based on the assumption cr^^^ = crf^j at energies in the range 546 < ^/s < 1800 GeV, we 
claim that this "enigma" will be solved by the forthcoming new very high energy colliders. The new 
Relativistic Heavy Ion Collider (RHIC) which has become operational at Brookhaven, should produce a 
value of (jf^j at y/s — 500 GeV very soon, allowing a comparison with the crfoj measured value at y/s — 546 
GeV. Later on, the LHC should give us a af^t value at y/s ~ 14 TeV. If it could be run at 2 TeV, a 
direct comparison between its a^^^ value and the a^^^ found at the Tevatron could be done. Additionally, 
experiments such as the HiRes |3 and the Auger Observatory |2] will bring new highlights to solve the 
puzzle about extrapolations to cosmic ray energies, i.e. whether the real values of atot are described by 
parameterizations consistent with a fast energy rise at high energies, which we may called "cosmic result" 
or by a slower rise in energy, the "accelerators result" as is the case presented in this paper. We argue 
that instead of looking for a new physics based on "exotic" events to explain very fast rise in energy of 
(Jtot, what we need is highly trustable data at intermediate energies (2 < y/s < 15 TeV) to use them 
as confident lever arms for extrapolations, to approach in this way to a more accurate description of 
atot at cosmic ray energies. In summary, extrapolations from accelerator data cannot be disregarded to 
constraint cosmic ray estimations. 

A CONFIDENCE INTERVALS OF THE MODEL PREDIC- 
TION 

The ability of a statistical method to reproduce data of any physical quantity with high precision gives 
the pattern for the prediction of out of the range data. In the context of pp and pp hadronic total cross 
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sections at very high energies a great deal of work has been done out of the energy range of accelerators: 
predictions are usually compared to cosmic ray data producing a disagreement which explanation has 
also been widely discussed in the literature. Such comparison depends critically of a highly confident 
band of uncertainty for any parametrization model. The validity of any statistical method to predict 
a given physical quantity out of the range values, "prognosis" (extrapolation) depends on its precision 
to reproduce the employed data (interpolation). A fundamental task of any prediction method is to 
minimize the error band of the predicted set of values. In the specific case of (Jtot what is searched is 
to obtain a prediction beyond the energy range of the employed data with the minimum of dispersion. 
In this context, popular techniques are derived from the statistical theory known as Regression Analysis 
either on the multiple regression approach or on the simplest version, a simple regression analysis, both 
based on the method of Least Squares. Among the important indicators of any statistical method within 
the frame of the present work there is the confidence level (C.L.) and the confidence intervals (e.g., 1491 ). 

For a given phenomenon characterized by an independent variable ( x ) that generates a response 
variable ( y) regression analysis supplies a procedure to estimate the corresponding statistics: to fit a set 
of data of the phenomenon in between the known points (interpolation), to estimate (jj) the mean value of 
(y) , to predict a future value of (y) beyond the known points {prognosis out of the data range) for a given 
value of (x) ), and to build a confidence interval around each of the estimated or prognosticated values 
(for instance |54l I55| ). In its general form the dependent variable (y) can be written as a function of k 
independent variables xi,X2, ■■■■Xk, where the variables Xj may represent powers of these variables, cross 
products of the variables, or even to present a parametric dependence of other variable (see for instance 
|56j). Statistic techniques of regression analysis are based on the minimization of the quadratic sum of 
data deviations with respect to the employed mathematical model of prediction: for a given distribution 
of residuals Ri = yi — yi a full set of confidence indicators are generated, where typical error bands (belts) 
for extrapolations up to m steps beyond the last n — esim experimental point is given as 

Error B = yn+m ± t^s/i^d (20) 

where 5^ is the variance, defined as the square of the standard deviation (usually called the standard error 
of estimate), yn+m is the corresponding Central prediction and ^1^/2'}^ denotes the Probability Density 
Fimction (p.d.f.) known as the Student's {t} -distribution for the k values of the independent variables 
with p degrees of freedom, and a confidence coefficient 5/2. However, the previous generalization does 
not incorporates effects related to the position of the experimental values of (y) around the employed 
Central Value, that is the correlation among residuals that generates the regression model |5()| . In other 
words, it is assumed no interaction terms, implying that each of the independent variables affects the 
response {y) independently of the other independent variables and the error e associated with any one {y) 
value is independent of the error associated with any other (y) value. Such an omission of the residuals 
correlation leads to a notorious modification of the statistic estimators [HHl- 

A.l The Forecasting Method 

The previous limitation, mentioned above, can be surmounted by identifying the obtained data set 

{xii,X2i,xzi Xki',yi) with a time series, and then to proceed an evaluation of the correlation among 

consecutive residuals (the so called autocorrelation procedure, from which the simplest one is the au- 
toregressive Ist-order model). The forecasting statistic method is based on autocorrelation models which 
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allows for the evaluation of the correlation among consecutive residuals, out and in the set of data by 
means of an iterative process |5()| . In general, this procedure modifies the fitting constants of the model, 
the estimated variance of residuals and minimizes the width of the error intervals for interpolation or 
prognosis (extrapolation), increasing in this way the confidence level {CL) of predictions j56| . To quantify 
the effect that has the autocorrelation of residuals on the regression model and associated estimators, let 
us represent such a model by a response variable y = E{y) + e where 

E{y) =Po + Pixi + [32X2 + + PkXk (21) 

represents the deterministic component of the proposed regression model, xi,X2, •■•,a;fe are the indepen- 
dent variables, that are known without error, and may represent higher order terms and even functions of 
variables as long as the function do not contain unknown parameters, /3i,/92, ■■■■,f3k are the unknown co- 
efficients to be determined by the least squares method, representing the contribution of the independent 
variables and e represents the random error component. For the evaluation of residuals Ri — yi — Eijji) 
it is assumed they have a normal distribution which mean is zero and the variance is constant. In the 
autoregressive model of Ist-order each residual Ri is related with the previous one as 

R, = (j)R,^i + r, (22) 

where (j> is the autocorrelation constant among the residuals (| (j) |< 1), |47l I48j . and in this case is 
a Residual called white noise, uncorrelated of any other Residual component. The incorporation of the 
effect of autocorrelation of residuals to the solution of the regression problem through equation (A2) 
leads to a modification of the regression constants and the corresponding variance: using the set of data, 
the estimation of the k + 1 regression constants of eq. (A2) and the constant cj) are obtained from the 
autocorrelation model according to the following interpolation equation for the response variable yf. 

Vi^ Pq+ PlXl^i + P2X2,i + PkXk,i + $Ri-i (23) 

where i — l...n, Xk.i represents the independent variable Xk corresponding to the point (i), the small 
hat indicates the estimated values, that is they have been evaluated, and the mean of the residual white 
noise has been taken as fn+i — 0. The prediction of values beyond the n-esim data value, beginning for 
instance with is then given as 

yn+l = /3o + /3lXi^„+i -I- /32a:2,r, + l + + l3kXk,n+l + Rn+1 (24) 

where Xk,n+i represents the independent variable Xk corresponding to the point {n + 1). The forecast 
(prognostic) for the response variable ijn+i is : 

Vn+l ^ Pq+ /3ia;i,„+i -I- i32X2.n+l + + PkXk,n+l + 4>Rn (25) 

Similarly, for a posterior value yn+2, 

yn+2 = l3o+ l3lXi^n+2 + /?2a;2,n+2 + + l3kXk.n+2 + (0) Rn (26) 

and so on successively, that is the forecasting of values out of the range of data is an iterative process 
where every new prognostic make use of the previous Residual. Unfortunately, as we move more and 
more far from the last data value, the potential error increases due to possible changes in the structure of 
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the regression model, or changes in the variance value as the prognostic are taking place. The later can 
be surmounted by estimating the variance associated with each prognostic (forecasting variance) through 
the correlation constant (jj; therefore, according to the autoregressive model of Ist-order within the data 
range up to the data n we have a constant variance S'j, and for one step out of the data range (i.e., n+ 1) 
the corresponding variance is S-j ^^^-^ — »S'j[l + (j)^], whereas for two steps beyond (n + 2) the forecasting 
variance is Sj = <S'y [1 + 0^ + 0"*] and for m steps beyond the range (n + m) the forecasting variance 
is = 5*1 [1 + (/)^ + 0^ + .... + (/)^(™^^)]. On this basis, for a prediction interval with a confidence of 

100(1 — S)% and a type t — student distribution, the amplitude of the prognostic intervals for (n + m) 
steps beyond the n — esim data value is given in 56 , as: 

Err or B - y„+„, ± t^J^ ^ S}[1 + dp"^ + (j)^ + .... + dp^^^^^-^)] (27) 

Due to the incorporation of the autocorrelation of residuals, the error bands evaluated in this way gives 
a higher Confidence Level than other methods of Regressive Analysis which ignore this effect HHj, what 
is translated in a width decrease of the prediction intervals as a consequence of the decrease of the 
estimated variance(e.g., Table 9.7 in jHE]). It must be emphasized that technically the meaning of the 
error bars in different methods derived from regression analysis is exactly the same , since they quantify 
the level of confidence (C.L.), that is, they express the probability that the "true answer" can be placed 
with 100(1 — 5)% of probability within the interval of aleatory measurements, or predictions, when other 
measurement is done under the same conditions. However, even if the concept is the same , every method 
predicts a different value of (C.L.): in the particular case of the forecasting, the essential point is that it 
uses the additional factor of the autocorrelation among residuals, what improves (1) the Central prognostic 
value and (2) the amplitude of the confidence intervals (error bars), whose narrowing is translated precisely 
in a high level of confidence '56'. 

A. 2 Matrix Approach in Regression Analysis 

This forecasting method is based on the multiple regression theory and consists in determining a prediction 
equation for a quantity y (dependent variable), that in turns depends on k independent variables {xi), 
that is 

k 

1=0 

(with fo{xo) = 1) , where fi are arbitrary functions of Xi, and j3i are the fitting constants. In the 
generalized version the variable Xi may depend on other parameters, i{s,t,..). Therefore, the 

application of a multiple regression model to a given problem leads to trace a system of n equations with 
n incognitos, so that its solution is better obtained through a matrix formalism. Denoting with Y the 
matrix of {n x l)-dimension of the dependent variables and with X the matrix of [nx {k + l)]-dimension 

of the k independent variables, the row l,xii,xi2, ,Xik multiplied by the column matrix of the (3's 

determines the value yi of the dependent variable, the row 1, 2^211 2;22, ,X2k multiplied by the column 

matrix of the /3's determines the value y2 and so on. 
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The variables contained in the matrixes X, Y can be related by the matrix equation Y — XB , which 
is the matrix expression of the prediction equation (A9). The [{k + 1) x l]-diniension matrix B contains 
the values of the constants j3i needed to write in explicit form the prediction equation (A9). The /3's 
constants can be determined by the least squares method 56 through the condition 



[yj — E{yj)Y = i?j = minimum 



(29) 



where yj is the j — esim measurement of the response variable and E{yj) is the estimated Central Value 

with Eq. (A9). The condition (AlO) is satisfied when = 0, (z = 1 k) which leads to a 

system of n equations with k{~ n) unknowns. This system in matrix formalism can be written (e.g., 
[nZI) as 

{X^X)B = X^Y (30) 

where X* denotes the transposed matrix of X and B is the matrix of the expected values of the (3' s. 
From (All) we obtain the solution of eq. (AlO): 



B= [x^xy^x^Y 



(31) 



where (A"*A")^^ denotes the inverse matrix of X^-X. Essentially this equation minimizes the quadratic 
sum of the deviations of points yj with respect to the prediction equation (A9) ( 56 p. 783). With the 
previous matrixes several statistical estimators are easily determined, such as the Sum of Square Errors 
(SSE) 

SSE ^Y*Y ~B\X^Y) (32) 
and the variance required to evaluate the confidence intervals is: 

SSE 



Sf = 



(33) 



where the denominator defines the number of degrees of freedom for errors, given by the number of /3i - 
parameters. Once the Central values are known, we then evaluate the confidence interval for a particular 
value of the response variable, 2/p, using the matrix of the particular values of the independent variables 
which determine the estimated value of y^. Such a matrix, namely A., denotes the column-matrix of 
(fc + 1) X 1 dimension, which elements {1, a;ip, X2p, . . . , a^fcp} correspond to the numerical values of the (3i 
appearing in equation (^49). Therefore, according to the associated statistic the confidence interval for 
prediction within the range of data is determined as (^^ P- 795): 



InterpB — y ^ ^ 



^8,2 ^SfA^iX^X) 
and for prognosis out of the data range as ([SS] p. 800): 



ExtrapB 



y ± h/2 



02 



AHx'xr^ A 



(34) 



(35) 
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Here y denotes the central prediction, A* is the transposed matrix of A. InterpB{+), ExtrapB{+) 
and InterpB{—), ExtrapB{—) denote the corresponding Upper and Lower bounds respectively. t\^s/2} 
denotes Student's {t}-distribution for the n values of the independent variables with p degrees of freedom. 
Estimation have been done with a precision of 100(1 — S)%, assuming 5/2 — 0.025, which corresponds to 
a value of 95%. 



References 



[1] 

[2] 

[3] 
[4] 
[5] 
[6] 

[7] 
[8] 
[9 
[lO: 

[11 

[12 
[13 
[14 
[15 

[16 
[17 
[18 
[19 

[2o; 



Engel, R., Gaisser, T.K., Lipari, P., Stanev, T., Phys. Rev. D 58, 014019 (1998). 
See http://sunshine.chpc.utah.edu/research/cosmic/hires/| 



The Pierre Auger Project Design Report. Fermilab report (Feb. 1997). 

Lipari, P., e-print Los Alamos National Laboratory arXiv:hep-ph/0301196 vl (12 Jan 2003). 



Donnachie, A and Landshoff, P.V., Phys. Lett. B 296, 227 (1992). 

Augier, C. et al., Phys. Lett. B 315, 503 (1993). 

Martini, A.F. and Mcnon, M.J., Phys. Rev. D 56, 4338 (1997). 

Amaldi, U. et al., Phys. Lett. B 44, 11 (1973). 

Amendolia, S.R. et al., Phys. Lett. B 44, 119 (1973). 

Bozzo, M. et al., Phys. Lett. B 147, 392 (1984). 

Amos, N. A., Avila, C, Baker, W.F., Bertani, M., Block, M.M., Dimitroyannis, D.A., Phys. Rev. 
Lett. 63, 2784 (1989). 

Abe et al., Phys. Rev. D 50, 5550 (1994). 

Avila, C. et al., Phys Lett. B445, 419 (1999). 

Matthiae, G., Rep. Prog. Phys. 57, 743 (1994). 

Proc. International Conference and 8th Blois Workshop on Elastic and Diffractive Scattering (EDS 
99), Protvino, Russia (27 Jun - 2 Jul 1999). 

Froissart, M., Phys. Rev. 123, 1053 (1961); Martin, A., Nuovo Cimento 42, 930 (1966). 

Augier, C. et al., Phys. Lett. B 316, 448 (1993). 

Lukaszuk, L. and Nicolescu, B., Nuovo Cimento Lett. 8, 405 (1973). 

Gauron, P., Lipatov, L. and Nicolescu, B. , Phys. Lett. B304, 334 (1993). 

Amaldi, U. et al., Phys. Lett. B 66, 390 (1977). 



23 



[21] Honda, M., Nagano, M., Tonwar, S., Kasahara, K., Hara, T., Hayashida, N., Matsubara, Y., Teshima, 
M., Yoshida, S., Phys. Rev. Lett. 70, 525 (1993). 

[22] Baltrusatis, R.M., Cassiday, G.L., Elbert, J.W., Gerhardy, P.R., Ko, S., Loh, E.G., Mizumoto, Y., 
Sokolsky, P., Steck, D., Phys. Rev. Lett. 52, 1380 (1984). 

[23] Baltrusatis, R.M. et al Proc. 19Th ICRC, La JoUa, 7, 155 (1985). 

[24] Luna, E.G.S. and Menon, M.J. On the total cross section extrapolations to cosmic-ray energies, 
e-print of The Alamos National Laboratory arXiv: |hep-ph /0105076| vl (2001). 

[25] G. L. Pryke, A Comparative Study of the Depth of Maximum of Simulated Air Shower Longitudinal 
Profile, ,astro- ph/0003442 (2000). 

[26] A. M. HiUas, Nucl. Phys. B (Proc. Supp.) 52B, 29 (1997). 

[27] R. S. Fletcher et al., Phys. Rev. D 50, 5710 (1994). 

[28] Glauber, R.J., Lectures in Theoretical Physics, Reading: Interscience, N.Y. (1956). 
[29] Glauber R.J., Matthiae G., Nucl. Phys. B 21, 135 (1970). 
[30] Nikolaev, N.N., Phys. Rev. D 48, R1904 (1993). 

[31] Gaisser, T.K., Sukhatme, U.P., Yodh, G.B., Phys. Rev. D 36, 1350 (1987). 
[32] Glauber, R.J., Velasco, J. Phys. Lett. B 147, 380 (1984). 

[33] Velasco, J., Perez-Peraza, J., Gallegos-Cruz, A., Alvarez-Madrigal, M., Faus-Golfe, A. and Sanchez- 
Hertz, A. Proc. 26^^ ICRC, Utah, 1, 198 (1999). 

[34] Perez-Peraza, J., Gallcgos-Cruz, A., Velasco, J., Sanchez- Hertz, A., Faus-Golfe, A. and Alvarez- 
Madrigal, M., International Workshop on Observing ultra high energy cosmic rays from space and 
earth. Proceedings of the AIP No. 566, 343 (august 2000). 

[35] Block, M.M. and Gahn, R.N., Rev. Mo. Phys.57, 563 (1985). 

[36] Bourrely, C., Soffer, J. and Wu, T.T. , Nucl. Phys. 247B, 15 (1984); Phys. Rev. Lett. 54, 757 (1985). 
[37] Garboni, G. et al., Nucl. Phys. B 254, 697 (1985). 

[38] Gudell, J.R., Ezhela, V., Kang, K., Lugovsky, S. and Tkachenko, N., e- pint from Los Alamos National 
Laboratory arXiv: hep-ph/9908218| V (13 Oct. 1999). 

[39] Menon, M.J., Can. J. Phys. 74, 594 (1996). 

[40] Rubinstein, R. et al., Phys. Rev D 30, 1413 (1984). 

[41] Schubert, K.R. in Tables on Nucleon-Nucleon Scattering, of the Landollt-Brnstein New Series of 
Numerical Data and Functional Relationships in Science and Technology, Vol 1/ 9a, ed. by K.H. 
Hellwege, Springer Verlag, Berlin (1979). 



24 



[42] Bozzo, M. et al., (UA4 Collaboration), Phys. Lett. 147B, 385 (1984). 
[43] Bozzo, M. et al., (UA4 Collaboration), Phys. Lett. 155B, 197 (1984). 
[44] Amos, N., et al., (E710 Collaboration), Phys. Lett. 247B, 127 (1990). 

[45] Cheng, H. and Wu, T.T., Expanding Protons: Scattering at High Energies, MIT Press, Cambridge, 
MA (1987). 

[46] Henzi, R. and Valin, P. , Phys. Lett. 132B, 443 (1983); 160B, 167 (1985). 
[47] Fuller, W.A. Introduction to Statistical Time Series, Wiley, New York (1976). 
[48] Chatfield, Ch. The Analysis of Time Series, Chapman and Hall (1996). 

[49] Montanet, L. Barnett, R.M., Groom, D.E., Hikasa, K. et al., Phys. Rev. D 50-3 1271 (1994). 

[50] M. Albrow et al.. Comparison of the Total Cross Sections Measurements of CDF and E811, Fermilab- 
TM-2071 february 1999. 

[51] Bueno, A. and Velasco, J. Phys. Lett. B 380, 184 (1996). 

[52] Avila, R.F., Luna, E.G.S. and Menon, M.J. , e-print Los Alamos National Laboratory 
|ajXiv:hep- ph/0105065 vl (7 May 2001). 

[53] Caravalho, P.A.S., Martini, A.F. and Menon, M.J., e-print Los Alamos National Laboratory 
|arXiv:hep-ph/0312243| vl (17 Dec 2003). 

[54] Cramer, H. Mathematical Method of statistics, Princeton Univ. Press, New Jersey (1958). 

[55] Eadie, W.T., Drijard, D., James, D., James, F.E., Ross, M. and Sadoulet, B. Statistical Methods in 
Experimental Physics, North Holand, Amsterdan-London (1971). 

[56] Mendenhall, W. and Sincich, T., A Second Course in Statistics: Regression Analysis, Prentice Hall, 
New Jersey (1993). 

[57] Basilevsky, A. et al., in Applied Matrix Algebra in Statistical Sciences, North Holland, New York 
(1983). 



25 



